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Abstract 

We show that the vacuum ground state energy for massive scalars on a 1-dim L-sites peri- 
odic lattice can be interpreted as the thermodynamic free energy of particles at temperature 1/L 
governed by the Arutyunov-Frolov mirror Hamiltonian. Although the obligatory sum over zero- 
point-frequencies is finite on the lattice, a renormalization prescription is necessary in order to 
obtain a physical sensible result for the lattice Casimir energy. The coefficients of every term 
in the large L expansion of the lattice Casimir energy are provided in terms of modified Bessel 
functions. 
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I. INTRODUCTION 



There are at least two main motivations for introducing a lattice structure in quantum 
models. In solid state physics one introduce lattice models as the Heisenberg spin chain 
or the Ising model to study the properties of matter built out of a granular atomic lattice. 
In the realm of particle physics lattice models based on continuum quantum field theories 
such as gauge theories, especially QCD, serve as an UV cutoff and they enable to study 
non-perturbative effects. 

The introduction of a periodic lattice with a finite number of L sites will affect the 
physical observables such as masses of particles or the vacuum ground state. Although the 
literature refers to influences of finite lattices as finite-size effects (see e.g. [ij), one should 
be aware that compared to the original QFT on continuous flat Minkowski space two steps 
of idealization are involved. 

The first one is compactification of space on a circle (or hypertorus in the case of several 
spatial dimensions) of size R. The corresponding finite-size effects on masses of particles were 
initially studied by Luscher [2j and the influence on the vacuum energy by imposing periodic 
boundary conditions 4>(x + R) = 4>(x) is sometimes called topological Casimir effect 0, 0]. 
Originally the Casimir effect [5] describes the phenomenon that a pair of perfectly conducting 
plates in the vacuum will feel a force attracting each other. This can be understood as a 
response of the quantum fields on the influence of imposed boundary conditions, which in 
the case of parallel plates are Dirichlet boundary conditions (For recent review, see e.g. [B]). 
In the case of periodic boundary conditions 0, U\, e.g. studying quantum fields on a circle 
or a cylinder, there is not literally a material boundary. In this sense we will use the terms 
'ground state energy' and 'Casimir energy' synonymously. 

The second step involves the discretization of the periodic space. One can expect that 
this will induce additional finite-sites effects on the physical observables. Actor et. al.jS] 
have investigated the Casimir effect on a lattice in the case of Dirichlet boundary condition. 
Their treatment of the lattice was as an alternative regularization procedure compared 
to mode number or momentum cutoff methods. Since their primary interest was in the 
continuum results, they ignored all 'lattice quantities' because they 'take care of themselves' 
by vanishing for L — > oo. A numerical study on the influence of lattice Casimir effects on 
the dynamics of kinks in a discretized sine-Gordon model was done by Speight [5]. 
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In the case of (1 + l)-dim relativistic integrable theories Zamolodchikov [10] suggested to 
use the thermodynamic Bethe ansatz (TBA) [TT] by interpreting the ground state energy of 
the theory on a circle of circumference R as the free energy of the same theory on a infinite 
line at temperature 1/R [TJ] • To apply this strategy on determining finite size-effects on the 
energy of strings in AdS§ x S 5 as advocated by Ambjorn et. al [13], Arutyunov and Frolov 
[H] have to introduce a mirror theory, since the gauge fixed worldsheet-theory has lost the 
explicit relativistic structure: 



For a review on the applications of TBA in the AdS/CFT context, see e.g. [151 116] - 

In this work we want to take the situation with a finite lattice structure and peri- 
odic boundary conditions for a scalar field at face value and ask, what kind of energy 
Ecas(L, R,m) can be attributed to the corresponding lattice vacuum state. Any regular- 
ization or renormalization procedures, if necessary on a lattice at all, should take place for 
fixed L. As we will argue in the course of this paper, we will indeed need a renormalization 
condition also on the lattice although all quantities by construction are mathematical finite 
at first sight. 

As ingredients we only use canonical quantization and zeta function regularization [T51 
[19] . adapted to the situation of a periodic lattice, to derive a compact analytic result for 
the Casimir energy. The expression can be interpreted as the free energy of particle at 
temperature 1/L governed by the mirror-typ Hamiltonian H m irr in ([l]), where the 'coupling' 
constant is given by the dimensionless combination of the physical parameters g = L/ (Rm). 
From this we can extract analytic results for the finite-sites corrections to the 'classic' Casimir 
energy E ~ 1/R in a large L expansion. Further, our result contains the leading Liischer-typ 
corrections E ~ e ~ mR anc l a 'wrapping'-type [13] correction E ~ g 2L as special cases. 

We tried to make the paper self-contained as an lattice approach to the Casimir energy, 
but think it is necessary to point out some formal similarities with results discussed in the 
AdS/CFT context. The plan of the paper is as follows. In Section 2 we review the general 
properties of lattice scalar fields and outline the shortcomings of a naive identification of 
the zero-point energy, while in Section 3 we will use a careful treatment by introducing 
the lattice zeta function associated to the lattice spectral problem to extract a physical 
sensible result for the lattice Casimir energy. In Section 4 we demonstrate that our result 




(1) 



3 



contains the well established continuum Casimir energies for scalar fields as limiting cases 
and a systematic large L expansion can be obtained. Finally in Section 5 we clarify some 
mathematical points which lie behind our physical result. 



II. SCALAR FIELD ON THE LATTICE 



For completeness we will review shortly the description of scalar fields on a spatial lattice 
in the Hamiltonian formulation [8J and then discuss the corresponding ground state energy. 

Consider a scalar field with mass m on a circle with circumference R = 2n, which is 
governed by the Klein-Gordon equation 



(9 t 2 -^)0(t,x)+m 2 0(t,x)=O, 



(j)(t,x + R) = <f>(t,x). 



(2) 



A straightforward discretization of this second order equation on a periodic L-sites lattice 
is given by 

1 



rj- 



!> m - 20, + + m 2 (j)i = 0, 



I = 0,...,L- 1, 



M) = M*)- (3) 



By taking the lattice spacing rj = R/L to zero 77 — ^- one gets back the continuous equation 
0- 



1 (vac) 



for all Z = L — 1 and 



1 (vac) 



+ <&z(£) will be treated as 



The vacuum field configuration is just given by y x 
the fluctuations around this background field 4>i(t) 

quantum fields in the canonical way. 

Setting $j(t) = e int $i gives an eigenvalue problem for the difference equation 



1 



7] z 



Z = 0,...,L-1, $l = $o- (4) 



This problem is equivalent to diagonalizing the L x L-matrix 

/ -2 + m 2 r] 2 1 ... 1 
1 -2 + mV 1 



D 



1 

rf 





V 1 



1 -2 + m 2 rf 1 
1 -2 + m V J 



(5) 
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and can be solved by elementary methods of difference equation calculus [17J, giving the 
fluctuation frequencies for this lattice fields as 



a 



• 2 f nn \ 9 

sin ^— j + m , 



n = 0,...,L- 1. 



The corresponding eigenmodes 



are normalized according to 



L-l 



Vl 



e L 



(6) 



(7) 



y2<Pn,l4>* n ',l = -<W 



L-l 



1 



(8) 



n=0 



The lattice analogs of the quantum field operator and its canonical conjugate momentum 
operator with II;/] = can now be expanded in terms of the eigenmodes as 

L-l 



„ iUntj. » i „fflni J.* At! 

e </>n,za n + e </> nji a;j 



•>(«) = E 4f [ 

n=0 V^'n 



(9) 



where the creation and annihilation operators fit and a n satisfy canonical commutation 
relations [a n , a\] = S„^ and [a n ,a n ,] = [a^a*,] = 0. 
The lattice Hamiltonian 

L-l 

V 



(10) 



1=0 



can then be written as 



L-l 



L-l 



< a n + ^ 



(11) 



n=0 n=0 

The L-sites lattice vacuum \0)l is defined by the condition that it is destroyed by all anni- 
hilation operators: 

d„|0) L = 0, n = 0, 1,..., L-l. (12) 



Acting with the Hamiltonian (11 ) on this state gives a bare value E for the vacuum energy 



H\0) L = E (L,m)\Q) L , 



(13) 
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with 

L-l 



E (L,m) = ^J2 n n, (14) 



n=0 



where £l n as in (|6|. Since the number of modes are finite and given by the number of sites 



L, there seems to be no ambiguity to consider the sum (14)as the vacuum or ground state 
energy of the lattice model: Ec as (L,m) = E (L,m). Mathematically, the finite number L 
of frequencies gives a finite sum with a finite result. There are two main arguments why 
this prescription cannot give the right physical answer. 

First of all, one expect that a lattice Casimir energy Ec a s{L, m) approaches for L — > oo in 
a smooth way the finite Casimir energy for periodic boundary conditions in the continuum 



as reported e.g. in (6J [7]. This cannot be satisfied by (14), since for L — > oo it becomes the 
notorious divergent series, which has to be renormalized. 1 

Further, Casimir energies should be considered as a special case of quantum corrections 
to the energy of classical background fields 0d[2T]. The background field in the case of the 
Casimir effect is just the vacuum field configuration C ; = 0( mc ) = with energy E c i(cf) c i) = 0. 
Only the fluctuations (fifi of the field 4> = 4> c i + 4>fi should be treated as quantum fields 
<f>fi — > <3>, regardless whether is a continuum field (f>(x) or a lattice field 0;. 

Now consider what happens with the frequencies ^ in the limit of large 'mass' m — > oo: 

Vt n ^m+— — sin 2 — (15) 



m 2n 2 \ L ) 

In the large mass limit even the frequency of the lowest lying fluctuation mode escapes to 
infinity. Since an infinite amount of energy is necessary to excite vacuum fluctuations it is 
physically reasonable, that no quantum fluctuations on top of the classical background field 
appear at all. Therefore one expect that also any quantum corrections to the energy vanish 
and one should be left only with the classical energy E = E c \ = 0. But a resummation of 



every term in ( 15 ) according to ( 14 ) gives 



2 "" ' m87r 2 
which implies E — > oo for m — > oo. 



1 I L 3 

E -)• -Lm+ —— - (16) 



In [8] the corresponding sum ( 14 1 for Dirichlct boundary conditions was treated by the authors as a lattice 



regularization of the divergent sum in the continuum and not as the physical lattice Casimir energy. 
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We see here very explicit that the demand for renormalization follows not from the 
appearance of a divergent infinite sum at the very beginning. Although in the discretized 
version of the massive scalar field only a finite number of modes contribute, an infinity shows 
up which needs to be renormalized 

E Cas = E ren -> 0, m -> oo. (17) 

In the case of free scalars this corresponds to a specific normal ordered Hamiltonian : H : = 
H — Esuik, where Esuik should contain all terms which have be subtracted. 



At least in the massless case m = one can resum the original expression ( 14 ) explicitly, 
which results in 

e ° = tM^l)- < 18 > 

Expanding this expression for L — > oo gives 

The first term is divergent in the continuum limit, the second term is the famous Casimir 
energy contribution, the third term is an example for the sub-leading contributions to the 
energy coming from the discretization. In the continuum limit we need a proper renormal- 
ization prescription to handle the divergent term proportional to L 2 and to end up with 

L ( 7T \ L 2 , . 

Bc "=S COt (2l)-^ (20) 

which gives for L — > oo 

ECas = ~ (21) 



III. SPECTRAL ZETA FUNCTION FOR SCALAR FIELDS ON THE LATTICE 

The previous discussion has shown that we will need a more careful analysis in order to 
determine the lattice Casimir energy. In the following we will write \ n = fl^ with Q n as 
in In order to guarantee that no zero modes are in the spectrum, we assume further 
m > 0. To proceed we introduce the spectral zeta function for the corresponding difference 
operator Q as 

Us) = » 1+2s J2K s (22) 

n=l 
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and the 'regularized' energy is given by 



E reg {s) = ^Cd(s), (23) 



where parameter /i is introduced to keep track of dimensions. The expression (23) is finite 
for Re(s) > in the limit m — > oo and for Re(s) > \ even in the limit L — > oo. We 
emphasize, that we don't need the regularization because the sum over lattice frequencies 
is a divergent series. This is not the case for the discrete field on a finite lattice. We need 
the regularization to handle the infinity in the large mass limit m — > oo. In the previous 
section we have evaluated the zero-point energy as sum over the explicitly known fluctuation 
frequencies. Now we will use an approach, which incorporate the fluctuation spectrum in a 



more intrinsic way, defined as the roots A (A) — 2 = of the spectral discriminant 



A(A) = 2 cosh 



L arcosh ^1 — JY^ ~ m2 )^ 



(24) 



Then we can rewrite the spectral zeta function in terms of a contour integral representation 
[TS] as 

CdOO = ^V +2s J dA A-^ ln(A(A) - 2) = ^/i 1+2s J dA A~^(A), (25) 
The resolvent 

fl (A) = |-.n(A(A)-2) = ^L (26) 

has therefore L poles on the positive real axis at the places of the eigenfrequencies and the 
integration contour 7 encircles them. The advantage of this representation is, that it can be 
applied also in cases where the fluctuation spectrum is not known. The only requirement is 



that one has a spectrum defining equation, such as the spectral discriminant in (24). 

Deforming the contour to lie along the branch cut on the negative real axis A G [—00, 0], 
one arrives at the following result: 

E reg (s) = E Bulk {s) + E Cas {s) (27) 



2 Obviously, the regularized energy is a function of zeta function argument s, the lattice sites number L, 
and the mass parameter m: E reg (s, L, m). In the following we will only mention the arguments explicitly, 
which are the most relevant ones in the course of arguing. 
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with 



EBulk(s) 
Ecas(s) 



Sin TVS 

L -u 

sin(7rs) 



, 00 r5 
l+2 - 1 dAA~ s — g(-A), 



7T 



l+2s 



d 



OA 



(28) 



and 



2tt 2 

g(A) = arcoshll -— (A-m 2 )J . (29) 

The splitting of the integral into two parts is valid for < Re(s) < 1. Since from the 
behaviour of the integrands one can estimate that the integral in -EWfc(s) is convergent for 
< Re(s) < 1 and the integral in Ec as (s) is convergent for Re(s) < 1. 

Inserting 

(30) 



<9 7T 1 

dX L ^/(A + m 2 )(l + g(A + m 2 )) 



into the integrand of Esuikis) and introducing the new integration variable x 2 = A + m 2 the 
integral can be solved for Re(s) > in terms of a hypergeometric function: 



Esuikis) 



■ -2S 1+2S7- 771 



' 2' ' 



ir 2 m 2 



(31) 



Although the integral in (28) was divergent for s < 0, we can use the result (31) as analytic 



continuation of Esuikis) to negative values of s and evaluate it also for s 



1. 
2 • 



E 



Bulk { 



■1/2) 



Lm, 



-E 



7T 



L 2 



ir 2 m 2 



(32) 



where E is the complete elliptic integral of second kind. Now we expand for L — > 00 and get 



E 



Bulk 



-1/2) 



L 2 „rm _L P (2) 1 

j2 ^Bulk 



with 



E 



(0) 
Bulk 



71 



m 

T 



+ e (0) + 



L 4 " 



(4) 
y .B«Zfc 



+ ... 



- + 2 In 2 + In - - In m 

2 7T 



(2) 
Bulk 



1 A 

it m 
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3 L 

2 In 2 — In h In m 

4 7T 



(33) 



(34) 

(35) 
(36) 



In general, the terms E^ k are 0(m 2n+2 ). From these expansions one can see, that E^uik 
encapsulates the terms, which diverge for m — > 00 and/or L — > 00. Applying the large mass 



subtraction scheme the renormalized lattice Hamiltonian (11) should be defined by 



:H:=H-E Bvlk (-l/2,L,m). 



(37) 



Acting on the vacuum state |0)l provides now 



:H : \0) L = E C as(L,m)\Q) L , (38) 



where Ec as (L,m) is given in (28). 



This renormalization ensures that for any given values L of lattice sites the ground state 
energy in the limit m — > oo vanishes and, as we will see in the next section, this choice is 
also compatible with the results known from the continuum theory. 

There are several integral representations for Ec a s{L,m) possible. After integration by 
parts one can write the Casimir energy of a free scalar field with mass m on a periodic lattice 
with L sites as 

1 f°° 

E Cas (L,m) = - / dKp'(K)\n(l-e- £ ^), (39) 



with 



P(k) = V k 2 — m 2 , e(k) = 2Larsinh • (40) 



A more illuminating representation can be obtained by a minor rescaling of the integration 
variable p = mp and introducing the dimensionless parameter g = Then one can write 
the Casimir energy more precisely as 

E Cas (L,g,m) = 2mf(L,g), (41) 

where the dimensionless scaling function 

/(4^) = ^ °|ln(l-e-^)) (42) 

can be interpreted as thermodynamic free energy of particles on an infinite line at a tem- 
perature \ governed by the Hamiltonian 

H g (p) = 2 arsinh QgVP 2 + M • (43) 

This rather exotic looking Hamiltonian is formally exactly the same as the one in the so 
called mirror theory, introduced by Arutyunov and Frolov by a double Wick rotation of the 
light-cone gauge- fixed AdS^ x S 5 worldsheet theory [H] . The emergence of this Hamiltonian 
related to a lattice model is not as surprising as it may seem at first sight. The eigenfrequen- 
cies ^ of the lattice field modes mimic the dispersion relation H w , s in obtained for the 
string worldsheet modes after light-cone gauge-fixing and are also similar to the dispersion 
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relation of giant magnons, solitonic states propagating on the worldsheet [22J. This simi- 
larity had previously led to musings about an underlying lattice structure of the AdS 5 x S 5 
worldsheet [23J. 



IV. DISCUSSION 

The spectrum of free scalars with bare mass m on a L-sites periodic lattice with circum- 
ference R starts with (g = L/(mR)) 

• the vacuum state 

|0) L , E = E C as = 2mf(L,g), (44) 

• the one-particle states 



4|0)l, E = J^j sin 2 + m 2 + 2m/(L,s), ra = 0,...,L-l 

(45) 

By investigating certain limits in the dimensionless parameters L and g, we are able to 



show that our result (42) contains some well known as well as some new results as special 
cases. 



The continuums limit 



The limit L, g — >• oo, while Rm = L/g held fixed, corresponds to the continuum limit of 



a scalar field on a compact circle with circumference R. From (42) one obtains an inverse 
power-law behaviour in large L: 



E Cas (L,R,m) E<®(R,m) + ^E<${R,m) + ^E^(R,m) + 0{L- % ) 
where the first coefficients have the following integral representations: 



71 



E {2 \R,m) 
E {A \R,m) 



R 3 



24vr 
R 



d« 



K 



rn 

5 roo 



Vk 2 - m 2 (e RK - 1) 



5760tt 



27 



2Rne RK 



I) 2 



(46) 



(47) 
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' * T • » • . 



10 



20 



30 



*:::::::::::: i i i i i ■ T 



40 



50 



FIG. 1: Relative deviation of the lattice Casimir energy Ec a s(L,27r,m) compared to the contin- 
uum result E^°'(2tt, m), evaluated for m = 9/10, 1/2 and 1/10 (blue, magenta and yellow dots, 
respectively) as a function of lattice sites L. This shows that lattice effects are under control also 
for a small lattice as long as mR < 1. 
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FIG. 2: Casimir energy depending on mass parameter m for fixed lattice L = 30 

For the leading order contribution it is preferable to introduce in the integral rep- 
resentation a rapidity parameterization 3 p = sinh(6 l ). Then one can see that this term can 



There exist also the possibility to introduce an elliptic parameterization of the integrand in ( 42 ) in terms 



of Jacobi elliptic functions. But it seems, this will not lead to an integration in terms of 'classic higher 
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FIG. 3: The sub- leading lattice correction E^ 2 \R,m) for fixed circumference R= ^, jq-tt and ir in 
magenta, blue and yellow. 

be interpreted as the free energy of an ideal relativistic Bose gas [12J at temperature 1/R 



m r°° 

E {0) {R,m) = - d^cosh(^)ln(l -e- mi?cosh(e) ) , 
71 Jo 



(48) 



as well as the Casimir energy of a massive scalar field on a circle of circumference R with 
periodic boundary conditions [6j [7] 



E^(R,m) = -— -KAmRn). 

tt ' J n 



it * — ' n 

n=l 



(49) 



Further, in the limit R — > oo it reproduce the leading Luscher corrections [21 121] : 



m) 



71 



d6 cosh(fl)e' 



-mR cosh(6') 



m 



TT 



Ki(mR). 



(50) 



The sub-leading term in (47) provides the first lattice correction of 0(L ) and can also be 



expressed in terms of modified Bessel functions: 
m 2 R 



E^ 2 \R,m) 



24vr 



3 V —K 2 (mRn) - 6mR V -KAmRn) + m 2 R 2 V KAmRn) 



n=l 



n=l 



n=l 



(51) 



transcendental functions' such as the modified Bessel function -Ky(z) in the case of the hyperbolic param- 
eterization. 
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£(4) 




m 



FIG. 4: The sub-sub-leading lattice correction E^(R,m) for fixed circumference R= ^, j^tt and 
7r in magenta, blue and yellow. 

The general term E^ n ' (2rr , m) in the large L expansion is composed of integrals of the 
type 



I(a, b, c, d, R, m) 



l / '2 1 \ — h 

CLK{K — m ) 2 K 



3 Rk\c 



(52) 



(e RK - l) d ' 

These integrals can be expressed in terms of modified Bessel functions and their properties 
are summarized in Appendix A. 



The sub-sub-leading lattice contributions is given in terms of (52) as 



E&{R,m) 



R 5 



5760tt 



[27J(-1, 6, 0, 1, R, m) - 2iU(-l, 7, 1, 2, R, m)] 



(53) 



By iterative application of the reduction formula (A5) the final expression in terms of Bessel 
function series will be of the following type: 

7 oo 

E {4 \R,m) = ^^c u {R,m)a u {n)K u {mRn), (54) 

f=3 71=1 

with a v (n) a rational function in n and prefactors c nu (R,m). 

In Fig.l we have examined the relative deviation of the lattice Casimir energy 
Ecas(L,2n,m) to the corresponding continuum result E^°'{27i, m), which shows that for 
mR < 1 the lattice corrections are under control also for small lattices L < 10. 

In Fig.3 and Fig. 4 the sub-leading lattice contributions E^ and E^ are plotted as 
function of the mass m for three different values of the compactification size R = rr/2, 7rr/10 
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and it. Although the values are roughly of the same order of magnitude, one should keep in 
mind that these contribution are suppressed by 1/L 2 and 1/L 4 , respectively. 



B. The massless limit m — )■ 



The limit g — > oo and L held fixed corresponds to the limit of massless scalars on the 
lattice. We can immediately obtain the corresponding Casimir energy by setting m = in 



(39): 



Ecas(L) 



TT 



d/dn(l-e- 2Larsin Hr' 



By construction, this should be equal to the previous mentioned result (20) 



E Cas (L) = ^cot 



TT 2 



(55) 



(56) 



with continuums limit L — > oo expansion as 



Ecas — > 



1 



TT 1 



12 720 L 2 



+ 0(L 4 



(57) 



On the other hand we can evaluate the massless limit from the expressions (49) and (51 ) 
which were obtained by taking first the limit L — > oo: 



E (o) 

ren ' 



oo 

— T- 

ttR ^— ' n 2 

n=l 



Ei 2 i 



m 2 R 
' 24tt 



6 



ttR 
48 



C(2) 



TT 

6R' 



+ 



48 



R 2 m 2 m 2 R 2 m 2 R 2 



Y- 



J n=l 



AttR 



C(4) 



TT 



360i?' 



(58) 
(59) 



which is in perfect agreement with (57) for R = 2tt. Notice also the appearance of Riemann 
zeta function ((s). 



Although (39) was derived under the assumption m > in order to avoid zero modes 



in the fluctuation spectrum, we can treat m now just as an infrared regulator which we 



removed m — > in the final step towards (55). Further we recognized that for an numerical 



evaluation of (55) it is advantageous to use the logarithmic representation of the inverse 



hyperbolic function: 



Ecas(L) = i jT dKln M - (kk/L + yJ(™/Lf 



(60) 
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C. Decompactified lattice limit 

There are two different limit procedures, which result in a decompactification R — > oo. 

• Consider the limit L — > oo while g held fixed. This corresponds to R, L — > oo while 
7] = R/L held fixed, which is the situation of an infinite lattice with constant lattice 
spacing rj. The Casimir energy is 



Clearly, in this situation the lattice corrections are exponentially suppressed 0(e L ) 
and vanish completely in the strict L — > oo limit. 

• The other decompactification limit can be obtained by taking g — » and L held fixed, 
which corresponds to a 'dilute lattice' with infinite large lattice spacing rj — » oo. The 
Casimir function can be approximated in this case as 



Therefore the first contributions of the Casimir energy on an dilute L-sites lattice are 
delayed up to order g 2L . Interestingly, this coincide with the behaviour of wrapping 
contributions to anomalous dimensions of certain N = 4 SYM gauge operators, which 
are expected not to be captured by standard spin chain Bethe ansaetze [131 121] • 

D. The roadmap 

The different limit procedures can be summarized in the following roadmap, where we 
also sketch the qualitative behaviour of the vacuum energy in the major branches. Starting 
point is a massive scalar field on a periodic lattice of any size R with an arbitrary number 
of sites L: 




(61) 




(62) 
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continuum on circle 

E ~ EW{R) + L- 2 EW{R) 





finite periodic lattice 
eq. (45), (46) 



Casimir 





E 



R 



Liischer 



E ~ e 



-mR 




lattice on infinite line 





Wrapping' 



2L 



V. COMMENTS 



The observation that the lattice Casimir energy has a non-zero value at all can be traced 
back mathematically to the following issue. For s = —1, —2, —3, ... the lattice spectral zeta 
function 



L-l 



( D (s,L,m) = 



n=0 



L\ . 2 /7rn\ 2 
- sm ( T )+m 



m 



L-l -a 

EE 



n=0 fc=0 



k 



L \ 



2k 



1 sin 2fc m 



V L / 
(63) 

has a finite binomial expansion and there is no harm to interchange the two finite summa- 
tions. Using 

'2k\ L 



L-l 



E 



n=0 



snr' 1 ' f ~j'^J 



fceN, 



(64) 



k J 2 2fc ' 

gives then the following alternative representation of the lattice spectral zeta function in 
terms of a Gaussian hypergeometric function 



Cd(s, L, m) = Lm 2s 2 F 1 ( s, -, 1; 



L 2 



27—5 9 9 / ' 



-1, 



-2, -3,... 



(65) 



which is essentially a polynomial of order — s in the argument — ^ 2 . That the result should 
be a polynomial can also be seen from the fact that for s being a negative integer the zeta 
function is just the trace of the corresponding matrix D given in ^ and taken to the power 
— s: 

( D (s,L,m) = Tr(D- s ). (66) 
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Nevertheless, as the discussion of the previous section has shown, for general Re(s) < 1 



and s 7^ —1, —2, ... the identity (65) does not hold and should be corrected as 

L 2 



Cd(s, L, m) = Lm ' Zs 2 F 1 ( s, -, 1; 



2' TT 2 m 2 



+ S D (s,L) 



where the additional correction term for Re(s) < 1 has been found as 



(67) 



5 D (s,L,m) 



sin us 



71 



dx(x — m ) 



2\-s 



d_ 

dx 



, / , — Larcoshf 

In I 1 — e v 



(68) 



For general values of s it is not allowed to interchange the double sum in (63). If one 
performs this changing regardless better knowledge one has to correct the induced error 



by the additional term (68), which for s = — \ is essentially the Casimir contribution 
5 D {-l/2,L,m) = 2E Cas {L,m). 



VI. SUMMARY 



The aim of this paper was to derive an unambiguous result for the ground state energy 
aka Casimir energy of scalar fields on a L-sites periodic lattice. From the physics perspective 
we have pointed out that although all relevant quantities are finite on a lattice we need a 
renormalization prescription to extract a physical meaningful result for the lattice Casimir 
energy. We have argued that in order to have consistency with established results for the 
ground state energy in the continuum, one has to interpret the Casimir energy also on the 
finite lattice as quantum correction to the energy of a background field configuration. In the 
case of free fields this results in an L-dependent normal ordering prescription of the lattice 
Hamiltonian. 

The final result for the Casimir energy can also be read as the thermodynamic free energy 
of particles at temperature 1/L on a continuous line governed by the Arutyunov-Frolov 
mirror Hamiltonian [H]. This interpretation indicates that the considerations on ground 
state energies in (l+l)-dim relativistic field theories on a circle outlined by Zamolodchikov 
[TO] seems to be valid at least for free fields in a strict lattice description. Further, it is a 
strong support for our choice of finite renormalization of the lattice Hamiltonian, that in 
several limiting cases well established results such as Liischer corrections, massive Casimir 
energy or the so called 'wrapping' contributions are obtained automatically. 

4 For a similar issue in the case of double series representations of Epstein zeta functions, see [25] . 
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As mathematical origin of a non-zero lattice Casimir energy we have identified the issue of 
non-interchangeability of a double sum in the defining representation of the lattice spectral 
zeta function. The contour integral representation delivers an elegant way to determine this 
remainder term, cf. this result with an attempt to use an Euler-Maclaurin-type formula. A 
further advantage of the contour integral representation of the spectral zeta function used 
in this paper is, that it can in principle also be applied in situations where the fluctuation 
spectrum is not explicit known. It would be interesting to apply this strategy to lattice 
models with interacting fields, which allow for kink states [S]. 
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Appendix A: Bessel functions identities 

We consider integrals of the type 



I(a,b,c,d,R,m) = I d/t(/t 2 - m 2 )f k \ ^ 1 (Al) 



(e RK - l) d 

with parameters a, b, c and d chosen in order to have convergence. 
Some basic examples of these integrals includes 



I{-l,0,0,l,R,m) = ^KoimRn), (A2) 



71=1 

OO 



1(1,0,0,1,-8,771) = ^-y^-KAmRn). (A3) 
R z — ' n 

n=l 



In general, any convergent integral of type (Al) with b = or b = 1 is expressible as a 
series of Bessel functions: 



I(a,0,c,d,R,m) = T ( -I- ^ (-lW )^K^(mRn), 

(A4) 

I(a,l,c,d,R,m) = — r — 52(-l) n ( n , H j^K^mRn). 
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By partial integration one can show the following algebraic relation between different inte- 



grals (Af ) 



(a + 2) J(a, b, c, d, m) = — {b — 1)1 (a + 2, b — 2, c, d, m) — Rcl(a + 2, b — 1, c, d, m) + 

+RdI(a + 2,b-l,c+l,d+l,m). (A5) 



Therefore any (Al) with b > 1 can be reduced by a finite number of iterative applications 



of (A5) to a linear combination of Bessel function series given in (A4). 
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